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Abstract 

We find exact solutions for a universal set of quantum gates on a scalable can- 
didate for quantum computers, namely an array of two level systems. The gates 
are constructed by a combination of dynamical and geometrical (non-Abelian) 
. phases. Previously these gates have been constructed mostly on non-scalable 

systems and by numerical searches among the loops in the manifold of control 
parameters of the Hamiltonian. 
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1 Introduction 



Let H{t) = H(R(t)) be a time dependent Hamiltonian acting on an n dimensional 
Hilbert space, whose dependence on time is via some control parameters collectively 
denoted by R(t). Let this hamiltonian have a /c-dimensional degenerate subspace 
Vo(t) spanned by the instantaneous eigenvectors |2(t)), • • • \k(t)), with energy 

Eo(t). When the point R(t) moves around a loop in the space of control parameters, 
any state |^>(0)) £ Vq evolves into a state 

\m) ^u(t)\m) = e* m r\m)- w 

Here 

g i/3(t) ._ e -i/ * E (t')dt> ^ ^ 

is the dynamical phase where Eo(t) is the instantaneous eigenvalue of the subspace 
Vb and 

r := P(e$c A W- dR ) = T(e$c A V- dt ), (3) 

is the non-abelian geometric phase which is the holonomy operator associated with 
the anti-hermitian connection A given by 

A Kij (R) = (i(R)\JL\j(R)), (4) 

or 

A ij (t) = A„ ij (R)— = (i(t)\-\j(t)). (5) 

Note that the symbols P and T in the first and the second integral of © refer re- 
spectively to the path ordering and time ordering of the exponential around the loop 
C in the control manifold. The basic property of a general holonomy operator is that 
it is independent of the way the loop is traversed in the parameter space. In some 
special cases it depends only on a few basic geometrical properties of the loop, like 
its area. In our case which will be discussed in detail later, in which the connection is 
constant, the holonomy depends only on this constant connection (which is nothing 
but the tangent vector on the loop at t = 0), and the total time needed for traversing 
the loop. Note that in the above formulas, the parameter t does not necessarily point 
to time, although we use this word for explicitness. It points to any single parameter 
which parameterizes the loop C. 

In the special case when Vo is one dimensional, T is the abelian geometrical phase 
and identical to the well known Berry phase. 

This general scheme when applied to the field of quantum computation takes the 
name of holonomic or geometrical quantum computation and the unitary operators 
thus obtained are called holonomic quantum gates. The problem of exact holonomic 
implementation of quantum gates is of great interest in the field of quantum compu- 
tation 00EllHElEll3IEllHllini- This is due to the fact that holonomic quantum 
computation, being geometrical in nature has a degree of stability against a class of 
errors [21 El HI]- 
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In particular it is known that these gates depend only on the loop and not on the 
speed with which they are traversed. Moreover this stability is related to the robust- 
ness of such gates against small perturbations of the traversed loops and 
against various noises. 

In the past few years many theoretical proposals for holonomic implementation of 
quantum gates have been reported in the literature [1211131 17] and some of them have 
been realized experimentally ^2]- At present we can say that there have been only 
sporadic successes in overcoming one or the other of the many obstacles in the way 
of a successful implementation of holonomic gates. 

Among these problems the requirement of scalability is the most important one. Let 
us briefly discuss this issue. 

It is well known that any unitary gate can be constructed to arbitrary precision from 
a combination of a universal set of gates. There are many choices for this universal 

/ i i \ 

set. One choice is the set {H, P(4>), C(tt)}, where H = I J is the Hadar- 

mard gate, P{4>) = ^ o e # J IS ^ e P nase S a te and C(<p) = diag(l, 1, 1, e 1 ^) is the 

controlled phase-gate. If we can act by these universal set of gates on any single qubit 
or any two qubits of our scalable system, then we can say that holonomic quantum 
gates have been constructed on the scalable system. 

We stress that such gates should be implemented on a scalable system, which is a 
crucial requirement for any viable candidate for quantum computation. A scalable 
candidate of quantum computer takes the form of an array of two level systems, or 
qubits. The array of identical systems can of course have higher dimensions but only 
two of their states will play the role of our computational qubit. Since any operation 
between remote qubits can be divided into elementary logic gates on adjacent qubits, 
it is sufficient to enact the universal set of gates only on single qubits and two adjacent 
qubits. 

An essential property of a scalable system is that the two qubit gates be realized 
on the tensor product of the same space on which the single qubit gates have been 
realized. However most of the proposals of holonomic computation so far suggested, 
lack this property. For example it has been shown that by abelian holonomy or Berry 
phase, one can implement the one-qubit phase gate on a single spin subject to a time 
varying magnetic field, and the two-qubit conditional phase gate on a pair of of cou- 
pled spins [JJ. On the other hand the Hadamard gate which is needed to be added to 
the above set if we are to have a universal set of gates, requires non-abelian holonomy 
PP, and is much more difficult to realize than the other gates. In fact the proposal for 
its holonomic realization is based on a completely different system consisting of two 
degenerate qubits and two ancilla qubits O EESj ■ 

In other words the hadamard gate is not implemented on the same qubit on which 
the phase gate was implemented. 

Of particular interest to us are the models based on iso-spectral Hamiltonians 
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which is reviewed in section (J2J) . Although in this case one can easily calculate the time 
or path-ordered exponential, the determination of exact solutions (exact loops in the 
parameter space) which lead to a universal set of quantum gates is difficult. Such an 
approach has been followed in , but with a numerical search among the class of loops 
for finding the required loop for each member of the universal set. Moreover in the 
approach of [S] , the single qubit gates are constructed on two dimensional subspaces of 
a three dimensional space, (i.e. k = 2, n = 3) and the two qubit gates are constructed 
on four dimensional subspaces of a five dimensional space, (i.e. k = 4, n = 5). 
This construction has the drawback that a two-qubit gate is not constructed on the 
tensor product space of two qubits, a requirement which is highly desirable for scalable 
quantum computation. In a related work Niskanen , Nakahara and Solomaa |Hj employ 
a three state hamiltonian of the form 



to implement the phase and the Hadamard gate on the degenerate subspace Vq = 
Span{|0), |1)} and the controlled phase-gate on the degenerate subspace of the Hamil- 
tonian 



This construction can be generalized to iV— qubit case. The dimension of the full 
Hilbert space scales as 3^. The one and the two qubit hamiltonians of the proposal 
of [SI are rather abstract and we do not know of any concrete realizations in terms of 
spins or some other suitable observables. 

Remark: Perhaps it is not strictly correct to say that the proposal of |13j is not 
scalable. In fact in this proposal which implement the universal set {e' l8cry , P(4>), C(tt)}, 
one can achieve scalability by going through polynomially more steps and using an- 
cilla bits which do not destroy the scalability. The aim of this paper is to achieve 
scalability by encoding each qubit in the two lowest states of two adjacent spins in a 
spin chain as shown in figure (0). This is equivalent to using s N ancilla qubits for 
the N computational qubits. 

What we want to show in this paper is that one can indeed find exact solutions 
for universal gates on a scalable candidate for qubits. For this aim, we consider a 
spin chain, and encode each qubit into the Hilbert space of two adjacent spins. We 
then show that by moving around appropriate loops in families of is-spectral spin 
chains, (determined by the adjoint action of suitable operators) one can implement 
a universal set of one and two qubit gates on such a spin chain. Such gates are 
realized as a combination of dynamical and geometrical non-abelian phases, when the 
corresponding loops are traversed. For each gate there are a set of parameters in the 
form of axis of rotations and frequencies which when tuned suitably enact appropriate 
one and two qubit gates on single and two adjacent qubits. 




/ 
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Hq 2 qubit = jj l qubit ^ j + j ^ 1 qubit_ 



(7) 
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Figure 1: In this array of spins, the degenerate ground state of the two spins in a 
block, say spins 1 and 2 encode the two states of a qubit. Exact holonomic gates act 
on each qubit and adjacent qubits by suitable manipulations of spins. 

The structure of this paper is as follows: In section (J2J) we review briefly the 
holonomies derived from iso-spectral hamiltonians. In section (j3J) we describe our 
model in detail where we show how a universal set of gates can be implemented by 
a combination of dynamical and geometrical phases on the spin chain. We end the 
paper with a discussion in section l|4jl. 



2 Holonomies of iso-spectral Hamiltonians 

In view of the time ordering of the exponential a closed formula for the holonomy 
operator © can not be obtained for most connections and one has to resort to 
numerical methods for an approximate calculation of this operator. 
It would be much desirable if we could find connections whose holonomy could be 
calculated exactly. This obstacle can be partially overcome by restricting ourselves 
(and hence paying a price of a using a limited source of holonomies) to constant 
connections for which we have 



r = 




(8) 
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where T is the total time required for traversing the loop. 

In view of the remark following equation Ijljl. the parameter T is not necessarily the 
total time needed for traversing the loop, it is only the final value of the parameter t 
which parameterizes the loop, which in the sequel will be set to unity. 
Perhaps the simplest way for obtaining time-independent holonomy operators is to 
consider iso-spectral family of Hamiltonians. These are the hamiltonians which are 
of the form 

H(t) := e xt H e- xt , (9) 

where X is any anti-hermitian operator and H$ is the hamiltonian at time with the 
degenerate subspace Vq, spanned by the vectors |1), |2), • • • \k), i.e. 

H \i)=E \i), i = l---k. (10) 

From relation Q we find that at any time t, the instantaneous eigenstates will have 
a simple form 

\i(t)) =e xt \i), H(t)\i(t)} =E \i(t)), i = l---k. (11) 

In this case we will find from @ that A will be constant, namely 

(i\A\j) = (i\X\j), i = l---k. (12) 

It is important to note that X is an anti-hermitian operator defined on the full Hilbert 
space and A is an operator defined only on the degenerate subspace and this relation 
implies only that the projection of X on this subspace is equal to A, that is 

X\ Vo = A. (13) 

Therefore a large number of operators X can lead to the same connection. Mathe- 
matically different loops in the parameter space are specified by their tangent vector 
at the origin which is nothing but the operator X. 

By re-scaling the time variable so that T = 1, and taking into account the inevitable 
dynamical phase, we arrive at the final form of the operator U acting on the space of 
a single qubit, 

U = e- iEo e A . (14) 

3 Holonomic computation on a spin chain 

We take an array of qubits as shown in figure ^ so that the spins within each block, 
say the spins (1,2) interact according to the following Heisenberg hamiltonian: 

H = B{o\ z + <j 2z ) + Jot ■ 02, (15) 

and different blocks do not interact with each other or their interaction is so weak 
that we can consider them effectively non-interacting at time t = 0. Later on we will 
make these blocks interact in order to implement two qubit gates. It can be easily 
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verified that if we choose the magnetic field so that B = 2 J, then the ground state of 
the Hamiltonian will be doubly degenerate. In fact the spectrum of the hamiltonian 
is as follows, where the states \4>q) and |V>o) are the degenerate ground states and \4>x) 
and \(f) 2 ) are the first and the second excited states. 



!+,+>, E = 5J, 

7§(l+'-> + |- +»> E = J ' 

I-,-), £ = -3J, 



\fh) = ^(l+,->-|- +», E = -3J. (16) 

We take the code or computational qubits to be the degenerate ground states, namely 

io> ee |0 O > = 1-,-), ii) = m = ^(i+ -) - i- +»• (17) 

In spin notations where S and S z respectively denote the total and the z component 
of the two spins, the code qubits are |0) = \S = 1, S z = —1) and |1) = \S = 0, S z = 0). 
We suppose that these states are not hard to access and control experimentally. 
At low temperatures the two spins reside in the degenerate two dimensional subspace 
which is a desirable situation for the initialization of the computer. If these two states 
were among the excited states of the Hamiltonian rather than the ground state, we 
would have been faced with an extra problem of exciting the two spins to these states. 

Let the operator X be of the following form 

X = in- (wi^i +w 2 ^ 2 ), (18) 

which describes the rotation of the spins Si and S2 around the axis n with frequencies 

uj\ and u>2 respectively. 

From ffl we find the gauge potential to be: 



A),o = (0o|^|0o) = +u} 2 )n z , 

i 

i 

A 1A = (VoWo} = 0. 



= ((t>o\ x \i>o) = ~ U2)(n x + iny), 

Ai,o = (i/)o\X\</)o) = -^=(uji- io 2 ){n x - in y ), 



Therefore in this subspace the gauge potential will be given by the following op- 
erator 

A = i(r x a x + r y a y + r z a z + r z I), (19) 

where 

1 -1 -1 

r x = —{u 1 - uj 2 )n x , r y = — (ui - L0 2 )n y , r z = — (ux + u; 2 )n z . (20) 
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After the lapse of time T = 1 and acquiring the dynamical phase 3 J, the gate 

U = e i ( r ^ c 'x+ryCTy+r z cT z + {r z + , iJ)I) 

will act on the the space of single-qubit codes, |0) and |1) . In this form the gate 
v! is not general enough, since the overall phase it applies namely r z + 3J is not 
independent of the other parameters. However at the end of any loop we can stop 
changing the parameters of the Hamiltonian, and only pause for a time interval r. 
This lapse of time will add a phase 3 Jt to the above phase and we obtain a general 
unitary gate given by 

U = e *( r ^ CT ^+ r 'y CT y+ r z "z+( r z+ 3J X 1 + T )H) (22) 

In this way by combining dynamical and geometrical phases we can construct any 
single qubit gate on our code qubits, since the parameters r x ,r y ,r z and r are inde- 
pendent. It is only necessary to choose the parameters ux,U2,n and r appropriately. 
We should emphasize that the specific structure of the degenerate states, namely that 
it consists of a product state and an entangled state, has been vital in our ability to 
arrive at a general form of the holonomy A = X\y Q with a simple choice of the oper- 
ator X. Had these two states been product states, we should have used complicated 
and hence unjustified forms of the operator X to arrive at the same general result. 
Now let us explicitly construct the two single qubit gates in the universal set, namely 
the phase gate P{4>) and the Hadamard gate H. 



3.1 The phase gate 

The phase gate is defined as 

W)=(j % ) =e i i^ +i ^\ (23) 

where m is an arbitrary integer. Comparison with (|22|) shows that for this gate we 
should have 

r x = r y = 0, r* = -|. r z + 3J(1 + r) = | + 2rmr. (24) 

From 1)20(1 we find that the following choice of the parameters implements this 
gate: 

n = (0,0,1), j l = U j 2 = , t i 3J(i + r ) = fj, + 2mvr. (25) 

This is a rotation of both spins around the z axis with equal frequencies followed by 
a pause for a time interval r given as above. The freedom in choosing the integer m 
guaranties that the time lapse r can always be positive as it should be. 
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3.2 The Hadamard gate 

The Hadamard gate is 



H 



1/1 1 



y/2 V 1 - 1 J ' 
We now note that H can be rewritten as follows: 

H = -j=(a z + a x ) = : k • 7?, 



(26) 



(27) 



where k is a unit vector k = -^=(1,0,1). Since for any unit vector k, e 1 2 k = ik • a 5 ', 
we find that 

= -i ik • 7? = -ie'f k ^ = e-^e^ = e ^+^*+^)+^ (2g) 

Comparison with (|21jl shows that the Hadamard gate is produced when we choose 
the following parameters: 



7T 



TV = r. 



r y = 0, r z + 3 J(l + r) 



7T 



+ 2m7r. 



2 X 2 ■ -v- ■ 2 

Comparison with ()2U() determines the parameters of rotation as follows: 



(29) 



wi = |\/3, w 2 = 0, 3J(l+r) = 



2+l)+2m7r. 
(30) 



Thus we have constructed our single-qubit gates on our code space which contains 
the computational qubits. We now turn to the conditional phase gate to complete 
our universal set of gates. 



3.3 The conditional Phase gate 

The controlled phase gate has the following matrix form when the basis vectors of 
the two qubits are ordered as 1 0, 0), |0, 1), 1 1 , 0) and |1,1): 



( 1 



V 



exp 



e i<t> J 



\ 



\ 



+ 2mir) J 



(31) 



In view of the demanded scalability, this gate should act on the space Vq (g> Vq, 
where Vq is the space on which the single qubit gates act. The space Vq<8>Vo is spanned 
by the states of two qubits, namely 



: ; ? 1 



|0,0) = |0o>®|0o> 



|0,1) = |^ )®|^ ) = -^(|-,-,+, 



1 



) ? ) 



+)), 



s 



1 



|1,0) = |Vo> ® |0o) = -^11+.- » . 
|1,1) = |^o>® lV*>> = o (I +,-,+,-) 



+ 1- +,-, 



-1-,+,-,-)), 

|+, -,-,+) - |-,+,+, 



To implement this gate the operator X should act on the Hilbert space of two adjacent 
pairs of spins, namely the pairs (5*1,52) and (53,54). Combining the geometric and 
the dynamical phases we find that the gate which will be implemented on the two 
qubits is equal to 



jj2qubit _ e A+i(6J(l+r)) 



(33) 



where A = X\y Q ®y and we have used the fact that the energy of the degenerate 
subspace is now 6 J instead of 3 J. 
We now take the operator X to be of the form 



X = i(f> {a 2z CT3 Z + a 2z + 032 



(34) 



This operator couples the endpoint spins of the two neighboring blocks which hitherto 
were considered non-interacting. 
It is easy to verify that 

X\0,0) = 

X\0,1) = 

X\1,0) = 

X\l,l) = 

This will then lead to 



-uf>\0,0) 

-#|o,i> 

-#|1,0) 



(35) 
(36) 
(37) 
(38) 



A = X 



Vq®Vq 



I -icf) 



V 



-%4> 



-i(p 



(39) 



0/ 



In view of (|33[) if we now choose r so that 6J(1 + t) = <j) + 2mir, we will find 



A + i6J{l + r) 



V 



+ 2mi\) J 



(40) 



and hence the conditional phase gate C(4>) will be exactly implemented on the two 
qubits. 

This completes our derivation of exact holonomies for a universal set of gates on an 
array of qubits. 
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4 Discussion 



We have been able to implement a universal set of quantum gates on a scalable sys- 
tem, by combining appropriately the dynamical and non-abelian geometrical phases. 
Our system consists of array of two-spin blocks each of which is a four dimensional 
space with a two dimensional degenerate subspace encoding the computational qubits. 
With these universal set at hand, one can construct any other gate to a sufficient de- 
gree of accuracy. The crucial step in this direction has been an appropriate choice of 
a degenerate subspace of a physical system which should represent the computational 
qubits. This subspace should be so that the gauge connection projected on it by a 
simple operator X be general enough to represent an arbitrary general gate. The 
choice of hamiltonian, so that its degenerate subspace has one entangled state and 
one product state has been essential in this step. Moreover we stress that the con- 
ditional phase gate has been constructed on the tensor product of two such qubits. 
We should add that the experimental realization of such a proposal is a completely 
different problem and we do not claim that this proposal is superior to others as far as 
experimental realization is concerned. We only emphasize the exact and the scalable 
nature of the proposal and hope that following the basic idea of this paper, namely 
taking two-spin blocks for representing qubits, other researchers can proceed to more 
practical and experimentally viable proposals. 
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